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Abstract 

We calculate the strong coupling constants gp**p* n , where P** (D**,B**) is the 
1 + p-wave state, in the framework of the light cone QCD sum rules, and using 
these values of gp**p*Tr, we compute the hadronic decay widths for D** — > D* it and 
B** -> B* -it. 



*Permanent address: Institute of Physics, Azerbaijanian Academy of Sciences, Baku, Azerbaijan 



1 Introduction 



The main goal of the future c - r and B - meson factories is a deeper and more comprehensive 
investigation of the properties of the heavy mesons, containing charm and beauty quarks. 
In particular, search for the excited states of D and B mesons and their decay modes 
constitute one of the main research program of the above mentioned factories, and can play 
an essential role for understanding the dynamics of the excited states. 

In general, for the interpretation of the experimental data from heavy meson physics, we 
need to know the large distance, i.e., nonperturbative effects. For example, for the exclusive 
decays which can easily be measured experimentally, we need a more accurate estimation 
of the form factors and other hadronic matrix elements that are described by the long 
distance effects, thus one needs a method which takes into account the nonperturbative 
(long distance) effects. Among the different approaches in estimating the large distance 
effects, the QCD sum rule method H occupies a special place, since this method is based 
on the first principles of QCD and the fundamental QCD Lagrangian. 

In this work we use a version of the QCD sum rule, namely the light cone QCD sum 
rule. This method is based on the Wilson Expansion of the T-product of currents near 
light cone, in terms of different non-local operators. These operators are characterized by 
their twists rather than their dimensions. Matrix elements of the non-local operators in 
the variable external field are identified with a set of the wave functions of increasing twist, 
and replace the vacuum expectation value of local operators that appear in the traditional 
sum rule method. The form of wave functions are restricted by the conformal invariance of 
QCD. More about the details of this method and its applications can be found elsewhere 



in the literature M- 19 



In this article this method is used for the calculation of the strong coupling constants 
gp**P* n of the excited positive parity meson decays, where p**(p*) is the 1 + (1~) meson 
state. In sect. 2 we derive the sum rule for the gp**p* n coupling constants. Sect. 3 is devoted 
to the numerical analysis, where we also compute the widths of the 1 + meson decay to P* n 
state. 
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2 Light Cone Sum Rule for gp**p* n Coupling Constant 

According to the general strategy of QCD sum rule method, the coupling constant gp**p* n 
can be calculated by equating the representations of a suitable correlator in hadronic and 
quark-gluon languages. For this aim we consider the following correlator: 

= ij d 4 x e^ x (Tr(q)\T {d(x)^Q(x)Q(0) ll/l5 u(0)} |0> . (1) 

Here d'j^Q^Q'y^u) is the interpolating current for 1~ (1 + ) meson state, Q is a heavy quark 
(charm quark for D meson, and beauty quark for B meson case), p\ is the momentum of 
the 1~ meson. When the pion is on the mass shell q 2 = m^, the correlator function (1) 
depends on two invariants, p 2 and p\. In what follows we set m n = 0. 

First consider the physical (hadronic) representation of (1). Physical part of it can be 
expressed in terms of the contribution of the lowest lying resonances P** and P* in the 
corresponding channels 

IV = <«Q)P (Pi)\P (p)> pj _ mh p2 _ ^ • (2) 
The matrix elements entering in eq.(2) are defined in the standard manner: 

(P*\d^Q\0) = mp.fp.eM , 
(0\Q^ 5 u\P**) = m P ..fp..e^(p) , (3) 

where mp* (m P **), f P * (fp**) and e M (e^) are the mass, leptonic decay constant, and 
the polarization vector of the vector 1~ (1 + ) meson state. In general, the matrix element 
(itP*\P**) can be written as 

(n(q)P*( Pl )\P**(p)) = F ( e * e «) + F 1 (qe)(q^) . (4) 

In obtaining eq.(4) we have used the transversality condition, p±e = pe^ = 0, and p\ = p—q. 
The second term in eq.(4) gives negligible contribution to the decay width P** — > P*7r, in 
comparision to the first term, since it is proportional to 

/ A 4 \ 



2 2 
TYlp** 777 p* 
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where A = nip** — nip* . Therefore we shall neglect the second term in (4) and set 
gp**P* T = F . Using eqs.(2), (3), and (4) for the physical part, we get 

_ m P *mp**fp,fp**gp**p* n \ p^p v Px^Pxv {V\V)V\p,Vv \ 

^ (p 2 — mp**)(pi — mp,) 1 mp„ mp* m P *m P ** J 

At this point we would like to make the following comment. Since the vector current qj^Q 

is not conserved, it also couples to J p = + scalar mesons Po as well as the J p = 1~ vector 

mesons. Therefore the Po contribution should be taken into account in the sum rule and 

this addition introduces further uncertainities. In order to avoid the + meson contributions 

we must choose a structure that does not contain its contribution. 

Noting that the corresponding transition matrix element is given as 

{Q\<rfnQ\ p o) = fDom D op lfl , 

that is, only the structure ~ pi^ contains the + meson contribution, it follows from eq.(5) 
that we have only two structures, namely ~ g^ v and p^Pv, which do not contain the J p = + 
meson contribution. In our analysis we choose the structure ~ g^ v . We also perform 
calculations for the structure Pp,p v and the final results of both structures in predicting of 
gp**P* n are practically the same. 

After performing double Borel transformation over variables — p 2 and — p\ (see eq.(5)) 
for the physical part of the sum rule for gp**p n , for the structure of g^ v , we get : 

UPhVS = M^9p~p>*fp~fp>m P ..mp.e\ ^~ ) . (6) 

Now we turn out attention to the theoretical part of (1). In this calculation we will use the 
notation of the work ||16|| . After a lengthy calculation we get (after double Borel transfor- 
mation over — p 2 and — p?) 



Jl1 Ml Ml 



rn 



2 



m u + m d 



m Q M <Pp(u q ) 



1 



+ 2m 2 Q g 2 (u ) - -MVIK) + 2(M 2 + m 2 Q )(g[(u ) + G' 2 (u )) 
1 A^\ n r , f 1 ^ . ^|(«i,l-«i-a 3 ,a 3 ) 



1 „r r u fL-oa 

+ -M 2 / dax / da 

2 '-JO Jun—ai 



3 ~ 

u JUQ—ct\ 

u o da 3 



>UQ—ai a 3 
/ ua 3 r 1 

+ / |^||(^o - a 3 , 1 - u Q , a 3 ) - <p\\(u - a 3 , 1 - u Q ,a 3 )j 



a 3 
1 da- 



\(p\\(u , 1 - u , a 3 ) + <p\\(u , 1-Uo- a 3 ,a 3 )\ 
io a 3 L J 

+ (continuum contribution) . (7) 



The pion wave functions (p n (u), ipp(u), g%(u) and G 2 («), have the twists r = 2, r = 3, 
t = 4 and r = 4 respectively, and they appear in the matrix elements of nonlocal quark 
operators as shown below (see |J and 



(7r(g)|J 7M 7 5 «(0)|0) 



X 



^) f 1 due^ ux g 2 {u) , 



qx I Jo 



(8) 



(7r(g)|di7 B u(0)|0) = f du e iqux ip P (u) 



m u + m d Jo 



(9) 



G 2 (u) 



g 2 (v)dv . 



(10) 



The functions <f>{ a i) an d ^(^i) are the twist-4 wave functions and are defined in the following 
way: 



(ir(q) \d{x)^^ b g s G a i}{ux)u{{)) |0) 



9m 



qx / 



Qa gfa 



qx _ 



P«iVi(aj)e' 



igx (01+1103) 



and 



(7r(g) \d(x)^g s G a/3 (ux)u(0) |0) 



/ 2^ \ / X/j^ 

9/3 - ?q 0/J*T 



<y.r 



Pa^i(ai)e' 



iqa;(«i+tta3) 



qx 



(12) 



Vcti 



1 



and 



<iai<ia2<ia3<5(l — «i — «2 ~ ^3) • 



In (7) we set 



M 



M\Ml Ml , _ dcp 

Ml + Ml ' M ° ~ Ml + Ml V ~ ~d^ u=u ° 
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We omitted the path-ordered factor Pe^ l9a Jo du X ^ A ^( UX )}^ s j nce j n the Fock-Schwinger 
gauge x^Ap = 0, it is trivial. In numerical calculations, we choose the symmetric point 
Uq — \ which means that quark and antiquark have equal momenta inside the pion. At 
this point the subtraction of the continuum can be done by substituting 



e m 2 



_ 2 o 

e m 5 " — e m 1 



(13) 



at least for the twist-3 contribution |L6|. However, we use this substitution everywhere in 
(7) since higher twist contributions are negligible. Equating eqs.(6) and (7) and using (13) 
we finally obtain the following sum rule for the strong coupling constant gp**p-x : 



9p**Pwfp**fp* 



1 



mp**mp* 

2 



2 2 
m p„* + m p* 

2 M 2 



f n M e - e x 



TTh" 2 771 I 



m, 



2 f 1 ■ + jj^l(^(«o)+^(«o)) 



+ 



+ 



JuQ—ai 



u o da 3 



o a 3 
1 dots 
o a 3 



{<p\\(u - a 3 ,l - u ,a 3 ) - ip\\(uo - a 3 , 1 - w ,ct 3 )} 
{<f\\(uo, 1 - u ,a 3 ) + <p\\(u , 1 - u - a 3 ,a 3 )} . (14) 



From (14) it follows that, for the calculation of the value of the strong coupling constant 
gp**P* n , we need to know the leptonic decay constants of the P** and P*, fp**, and fp*. 
The decay constant fp** can be obtained from the two point sum rules: 



f2 2 

Jp**Trip,, 



57T J m 



•so 



dse 



(s - m\f 



2 + 



+ m Q (qq)e 



2 2 
m p**~ m Q 



AM 2 



(15) 



where 



(qq) 



(0.8 ±0.2) GeV 2 
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To obtain fp* it is necessary to make the following replacements: rap** — > mp* and 
change the sign in front of the second term in (15) (in this case, of course the value of 
the continuum threshold must also change ). Note that we do not take into account the 
perturbative 0(a s ) corrections in (15), as they are not included in (14) either. Note that 
the values of the decay constants f B * and fo* we use in our calculations are given in |16| . 



3 Numerical Analysis 



For the numerical analysis of the QCD sum rule (14) we first give the values of the input 
parameters: 

m 

f n = 132 MeV, - — (1 GeV) = 1.65 GeV, m c = 1.3 GeV, 

m u + m d 

m b = 4.7 GeV, m D ** = 2.420 GeV, m D * = 2.01 GeV, 
m B * = 5.279 GeV, m B ** = 5.732 GeV, 
(s ) D = 6 4- 8 GeV 2 , (s ) B = 35 4 40 GeV 2 . 

Using these parameters, from (15) for the leptonic decay constants /d** and f B ** we get 

f D „ = (300 ± 30) MeV , (16) 
/ B « = (200 ± 20) MeV . (17) 

Sum rule for gp**p* n contains the nonperturbative quantities, namely the wave functions. 



In our numerical analysis we use the wave functions proposed in (see also |]I6|). Having 
the values of the input parameters, one must find the region of Borel parameter M 2 , for 
which the sum rule eq.(14) is reliable. The lowest value of M 2 is usually fixed by imposing 
the condition that the terms proportional to the are resonably small. The upper bound 
for M 2 is usually fixed by the condition that the continuum and higher states contributions 
constitute about (25 4- 30%) of the ground resonance contribution. Under these conditions 
the fiducial range of M 2 for B (D) case turns out to be 8 GeV 2 < M 2 < 20 GeV 2 (2 GeV 2 < 
M 2 < 6 GeV 2 ). Using the values of the input parameters we get 

Ib** Jb* 9b** b*it = (0.78 ± 0.12) GeV 3 , (18) 
f D **fD*g D **D*ir = (0.68 ± 0.10) GeV 3 . (19) 
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Dividing this product by the decay constants, we finally obtain for the D**D*% (B**B*n) 
coupling constants : 

g B .. B .„ = 24 ± 3 GeV , (20) 
£Id**d*tv = 10 ± 2 GeV . (21) 

Substituting these values in the expressions for the decay widths, 



r(p 



o _^ p*+ _ 9p**p* n ^ _|_ ( m p** + nip*) 2 ^ ^ 



24tt 



Amp* 



{m 2 P ** — (mp* + m^) 2 } {nip** — (nip* — m n ) 2 } 



we get, 



T(D**° -> £>*+ 7T - ) ~ 249 Mel/ , 
r (_B**o _^ s *+ ^ 296 MeV . 



(22) 



(23) 



Strong coupling constants (and correspondingly the decay widths) of the decays P** — 
p *o p**- _^ p*- n o and p**o _^ p*o n o can eagily be obtained from P**° -> P* + tt" 

with the help of the isotopic invariance. 
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